Genetics and Molecular Biology, 33, 1, 51-56 (2010) 

Copyright © 2009, Sociedade Brasileira de Genetica. Printed in Brazil 

www.sbg.org.br 



Research Article 

Microsatellite analysis to estimate genetic relationships among five bulgarian 
sheep breeds 

Szilvia Kusza 1 , Doytcho Dimov 2 , Istvan Nagy 3 , Zsuzsanna Bosze 3 , Andras Javor 1 and Sandor Kukovics 4 

1 Institute of Animal Science, University of Debrecen Centre of Agriculture Sciences and Engineering, 
Debrecen, Hungary. 

2 Department of Animal Sciences, Agricultural University of Plovdiv, Plovdiv, Bulgaria. 
3 Agricultural Biotechnology Centre, Godollo, Hungary. 

4 Research Institute for Animal Breeding and Nutrition, Herceghalom, Hungary. 
Abstract 

Herein, genetic relationships among five breeds of Bulgarian sheep were estimated using microsatellite markers. 
The total number of alleles identified was 226 at the 1 6 loci examined. D A distance values were used for phylogenetic 
tree construction with the UPGMA algorithm. The two Tsigai and two Maritza populations were found to be 
geneticallvery closely related to each other y (0.198, and 0.258 respectively). The Pleven Black Head population was 
distinct from the other four. These results could be useful for preserving genes in these breeds, thereby ensuring 
their preservation in Bulgaria. 

Key words: genetic relationships, Bulgarian sheep breeds, microsatellite. 
Received: November 11, 2008; Accepted: October 10, 2009. 



Introduction 

Over the past decade, numerous studies on genetic di- 
versity in domestic livestock (mainly in small ruminants), 
based on the analysis of microsatellite loci, have been car- 
ried out worldwide. In earlier studies on sheep genetic 
structure, use was made of polymorphisms of phenotypic 
traits (coat color, wool color, horn type, tail length) or bio- 
chemical markers (blood groups, milk and blood protein 
variation) (Grigaliunaite et al. , 2003). The breeds examined 
were chosen mainly in view of their long history of isola- 
tion, unique phenotypic qualities or evolution within a 
singular environment. In old, isolated breeds, it can be as- 
sumed that uniqueness in ancestry and phenotype corre- 
spond. However, modern breeds with distinct, selected 
external characteristics may have become genetically simi- 
lar through gene flow that typically occurs through male- 
mediated crossbreeding, and the use of few commercial 
sheep breeds as the basis for breed development (Terrill 
and Maijala, 1991). Molecular genetic studies on popula- 
tion structure may improve the understanding of present- 
day genetic resources (Moritz, 1994). 

In Bulgaria, the differences among several breeds of 
sheep are related to the geographical location, viz., low- 
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lands, hills, or mountains, where they are raised, and there 
is evidence of a high level of genetic diversity (Kukovics 
and Javor, 2002; Dimov, 2006). In the lowland areas, dairy 
breeds, viz. , Pleven Blackhead, Patch-faced Maritza, White 
Maritza, Stara Zagora and various crossed populations, are 
the most frequent. In the mountains, Tsigai are typical, 
along with Karakachaska and a number of other breeds. 

The first Tsigai individuals appeared in Bulgaria be- 
fore the second world war. After the 1950's, Bulgarian 
mountain breeds were crossed with Soviet Tsigai variants, 
thereby giving rise to a new Bulgarian variant, consisting of 
two main types, the north-western (Staroplaninski Tsigai) 
and the south-western Bulgarian (Rodopski Tsigai). These 
breeds are completely white, thus differing from the origi- 
nal Tsigai populations which are multicolor or black to gray 
(Kukovics, 2006, Kukovics and Kume, 2006). In Bulgaria, 
the number of these Tsigai populations is decreasing. 
Kukovics and Javor (2002) stated that Rodopski Tsigai had 
8,000, while the Staroplaninski Tsigai had 32,000 breeding 
ewes in the country. However today, the population size of 
these two breeds is much smaller and they are considered as 
threatened with extinction (Dimov D, personal communi- 
cation). Milk production in Bulgarian Tsigai breeds is rela- 
tively low (50-60 L over a period of 100-120 days). This is 
a most likely consequence of earlier breeding being con- 
centrated on those traits involved in meat and wool produc- 
tion, then considered of greater importance. The other 
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breeds under study are used for milk production, with 
yields of around 100-120 L of milk over a period of 
120-150 days, with high priority also being given to lamb 
production. Some breeders and scientists have argued that 
Pleven Blackheads, predominantly bred in the northern re- 
gion of the country around Pleven, might not belong to the 
Tsigai group. As their milk yield is the highest, this makes 
them the most productive.. This breed also represents a sig- 
nificant part of the total breeding stock, with about 220,000 
individuals distributed among numerous farms in small 
flocks of 20 to 120 individuals (Dimov, 2006). 

The Patch Faced Maritza and the White Maritza 
sheep are separately bred in the region north of Plovdiv. 
These populations are small and could be considered under 
threat, consisting of only 1 ,000 heads per breed, with about 
500-600 ewes each, distributed on 14 and 15 farms, respec- 
tively. Projects targeting the preservation of these breeds 
commenced only recently, in the late 1990s (Dimov, 2006). 
Milk yield per ewe of Patch Faced Maritza and White 
Maritza is similar (1 12.76 and 101.60 L, respectively), and 
no significant difference is to be found in the prolificacy co- 
efficient (1.34 and 1.29, respectively) (Dimov and Kuz- 
manova, 2007). 

There are known differences in body measurements 
and weights, coloring, and even in performance among the 
examined breeds. The average body- weights and wither- 
heights of the breeds under study are summarized in Ta- 
ble 1. 

The aim of this study (which is part of a South, East 
and Central European extended project) was to investigate 
genetic similarities among different breeds of Bulgarian 
sheep, notably Tsigai, Pleven Blackhead and Maritza, in or- 
der to facilitate their rational development, utilization and 
conservation. 

Material and Methods 
Material 

Hair samples were taken from randomly selected in- 
dividuals from different herds (Table 4). 



In order to ensure that the animals chosen were unre- 
lated, the wool samples from White Maritza sheep were 
collected from five (n = 8, n = 11, n = 5, n = 12, n = 5), 
Patch-Faced Maritza from four (n = ll,n = 9, n= 10, n = 9), 
Pleven Blackhead from two (n = 19, n = 16), Rodopski 
Tsigai from two (n = 15, n — 15) and Staroplaninski Tsigai 
from two (n = 20, n = 22) different herds. The selected 
two-two flocks from Pleven Blackhead and the two Tsigai 
breeds were chosen from different parts of Bulgaria so as to 
avoid genetic similarities. 

The White Maritza and Patch-Faced Maritza flocks 
were at least 30 km apart. When pedigree information was 
available, those 3 to 5-year-old animals with no shared 
common ancestor for at least two generations were se- 
lected. Animals without pedigree information were ran- 
domly chosen, with the premise that all were 3 to 5 years 
old and known by the head of the Maritza Breeder Associa- 
tion, through wool samples certifying to minimal genetic 
similarities. Numbers of animals studied per breed are pre- 
sented in Table 4. 

Methods 

Genomic DNA was extracted from the hair samples 
as previously described (FAO/IAEA, 2004). Sixteen 
microsatellite markers (CSSM43, MAF35, TGLA53, 
TGLA357, INRA127, OarCP20, BM1314, OarAE119, 
MAF65, MAF70, MCM527, MCMA7, OarFCB20, 
ILSTS 1 1 , OarCP49, BM6506) were selected based on their 
level of polymorphism, location on different chromosomes 
and recommendation by one or more of the following orga- 
nizations: United States Department of Agriculture 
(USDA); Australian Gene Mapping Web Site; Food and 
Agricultural Organization (FAO); International Society for 
Animal Genetics (ISAG). PCR and genotyping protocols 
followed those of Kusza et al (2008). Arlequin ver.2.0. and 
Populations ver. 1.2.28 programs were used for data analy- 
sis (Weir and Cockerham, 1984; Langella, 1999; Schneider 
et al., 2000). The Nei (1987) minimum genetic distance 
(Da) was estimated to define the genetic differences be- 
tween populations. Phylogenetic trees were constructed us- 



Table 1 - Main differences in body size in the studied sheep breeds. 



Bulgarian sheep breed Adult rams Adult ewes 





Live weight (kg) 


Wither height (cm) 


Live weight (kg) 


Wither height (cm) 


Pleven Blackhead 


95-120 


74-80 


75-95 


68-78 


Staroplaninski Tsigai 


70-80 


65-75 


45-48 


58-65 


Rodopski Tsigai 


80-85 


65-75 


45-50 


55-65 


Patch Faced Maritza 


100-120 


75-82 


65-95 


70-80 


White Maritza 


100-120 


74-81 


70-95 


70-80 



Kukovics and Javor, 2002, Dimov D, unpublished results. 
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ing the Phylip Ver. 3.57c software package (Felsenstein, 
1995). 

Results 

Microsatellite loci 

All microsatellite loci were found to be polymorphic, 
with the number of alleles per locus ranging from 7 
(MAF35) to34 (MAF70). Genotyping of the MAF70 locus 
was repeated due to the high number of alleles, thereby con- 
firming correctness. In all the examined populations, the 
highest number of alleles was at this very locus. The total 
number of alleles was 226 at the 16 studied loci. The mean 
number of alleles per locus ranged from 4.8 (MAF35 and 
MAF65) to 15.6 (MAF70) (Table 3). Of all the 16 loci, 14 
population-specific loci were identified in the five breeds 
(Table 2). There were no population-specific alleles at the 
loci CSSM43, MAF35, TGLA53, TGLA357, INRA127, 
OarCP20, BM1314, OarAE119 and MAF65. Allele fre- 
quencies are available from the authors on request. 

Observed and expected heterozygozities are pre- 
sented in Table 3. At all loci, expected heterozygosity 
(H exp ) was higher than the observed heterozygosity (H obs ), 
ranging from 0.580 (MAF65) to 0.867 (MAF70), with a 
mean of 0.736 among the loc. H 0 b s ranged from 0.256 
(MAF65) to 0.713 (ILSTS1 1), with a mean of 0.523. 

Genetic structure of the populations was analyzed by 
Wright's F statistics. The F IS fixation index among loci var- 
ied from 0.051 (TGLA357) to 0.592 (OarAE119). The 
mean value of Fis was higher than Fst (0.288 > 0.083), and 
heterozygote deficiency was detected at all examined loci. 
The F ST suggested that around 8,3% of the total genetic 
variation was due to variation within studied populations 
and the remaining 91,7% corresponding to differences 
among individuals. 



Table 2 - Population specific alleles. 



Locus 


Allele 


Carrier individual from the population 


MCM527 


150 


Staroplaninski Tsigai 




159 


White Maritza 


MCMA7 


259 


White Maritza 


OarFCB20 


1 15 


Pleven Blackhead 


MAF70 


128 


Staroplaninski Tsigai 




173 


White Maritza 




131 


White Maritza 


ILSTS11 


296 


Staroplaninski Tsigai 


OarCP49 


107 


Staroplaninski Tsigai 




120 


Rodopski Tsigai 




124 


Pleven Blackhead 




89 


Pleven Blackhead 


BM6506 


184 


Patch Faced Maritza 




185 


White Maritza 



Table 3 - Mean number of alleles, observed (R abs ) and expected (H exp ) 
heterozygosities, and results of F statistics. 

Locus M* H obs H exp Fis Fst Fit 



M*: mean number of alleles per locus 

Population 

The mean number of alleles per population ranged 
from 6.25 (Rodopski Tsigai) to 8.625 (White Maritza) (Ta- 
ble 4). 

Mean observed and expected heterozygosities per 
population ranged from 0.458 (Staroplaninski Tsigai) to 
0.577 (Patch-Faced Maritza), and 0.726 (Staroplaninski 
Tsigai) to 0.798 (Pleven Blackhead), respectively. All pop- 
ulations had lower than expected levels of heterozygosity. 
Genetic diversity was highest in Patch-Faced Maritzas. Fis 
values ranged from 0.246 (Patch-Faced Maritza) to 0.376 
(Pleven Blackhead), and were positive in all the popula- 
tions. The heterozygote deficit was the highest in Pleven 
Blackhead and the lowest in Patch Faced Maritza among 
the examined populations (Table 4). 

From Nei minimum genetic distance (D A ) values, it 
was apparent that the two Tsigai (Staroplaninski and Ro- 
dopski) and the two Maritza (Patch-Faced and White) pop- 
ulations were genetically the closest among those studied 
(0.198, and 0.258 respectively) (Table 5). Pleven Black- 
heads were most closely related to Patch-Faced Maritzas 
(0.291) and most distantly related to Rodopski Tsigais 
(0.390). 

Figure 1 depicts the Neighbour- Joining dendogram 
constructed by way of D A values. Genetic distances among 
Bulgarian populations are apparent. The Tsigai and Maritza 
populations clustered separately in the tree, with the Pleven 
Blackhead more closely related to the Maritza sheep breeds 
than the Tsigais (closest to the Patch Faced Maritza). 



MAF35 


4.8 


0.475 


0.649 


0.268 


0.077 


0.324 


CSSM43 


9.4 


0.541 


0.781 


0.308 


0.046 


0.340 


MCM527 


7.6 


0.438 


0.772 


0.433 


0.068 


0.472 


TGLA53 


8.6 


0.550 


0.825 


0.334 


0.076 


0.385 


MCMA7 


8.8 


0.582 


0.778 


0.252 


0.098 


0.325 


OarFCB20 


6.6 


0.441 


0.732 


0.398 


0.084 


0.421 


TGLA357 


7.0 


0.665 


0.701 


0.051 


0.034 


0.083 


INRA127 


5.0 


0.625 


0.666 


0.062 


0.161 


0.212 


MAF70 


15.6 


0.440 


0.867 


0.493 


0.067 


0.528 


MAF65 


4.8 


0.256 


0.580 


0.559 


0.092 


0.560 


ILSTS11 


7.4 


0.713 


0.760 


0.062 


0.065 


0.122 


OarCP20 


5.2 


0.647 


0.702 


0.078 


0.139 


0.206 


OarCP49 


8.2 


0.630 


0.745 


0.154 


0.126 


0.261 


BM1314 


8.6 


0.533 


0.821 


0.350 


0.091 


0.410 


BM6506 


7.0 


0.536 


0.686 


0.219 


0.062 


0.267 


OarAE119 


5.0 


0.291 


0.714 


0.592 


0.088 


0.628 


Mean 




0.523 


0.736 


0.288 


0.083 


0.349 
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Table 4 - Observed and expected heterozygosities and Fis values in Bulgarian sheep populations. 



Population 


Samples per population 


Mean number of alleles per population 


Hobs 


H e xp 


Fis 


Staroplaninski Tsigai 


42 


7.2 


0.458 


0.726 


0.366 


Rodopski Tsigai 


30 


6.3 


0.520 


0.735 


0.290 


Pleven Blackhead 


35 


7.6 


0.490 


0.798 


0.376 


Patch Faced Maritza 


39 


7.8 


0.577 


0.769 


0.246 


White Maritza 


41 


8.6 


0.568 


0.787 


0.275 



Table 5 - Nei genetic distances (below diagonal) between the studied Bulgarian sheep populations. 



Staroplaninski Tsigai 


Rodopski Tsigai 


Pleven Blackhead 


Patch Faced Maritza 


White Maritza 


Staroplaninski Tsigai 


0.000 










Rodopski Tsigai 


0.198 


0.000 








Pleven Blackhead 


0.339 


0.390 


0.000 






Patch Faced Maritza 


0.490 


0.473 


0.291 


0.000 




White Maritza 


0.501 


0.619 


0.313 


0.258 


0.000 



Staroplaninski Tsigai 



' Rodopski Tsigai 

I Pleven Blackhead 

I Patch Faced Maritza 

White Maritza 

Figure 1 - UPGMA tree constructed from D A distances, showing the rela- 
tionships among the five Bulgarian populations. 

Discussion 

Recently, the preservation of unique, genetically dis- 
tinct breeds of domesticated animals, especially indige- 
nous, has received much attention. Knowledge and 
information on genetic diversity, population structure and 
genetic relationships between populations are absolute pre- 
requisites for defining and accomplishing effective preser- 
vation strategies. This study aimed to characterize the 
genetic diversity and structure of two Bulgarian Tsigai, two 
Maritza and the Pleven Blackhead sheep populations by us- 
ing sixteen microsatellites. 

For the sampling process, unrelated individuals were 
selected from various flocks, from different parts of Bul- 
garia, to so avoid genetic similarities among the animals to 
be studied. In all, 226 alleles were detected at the 1 6 loci un- 
der consideration. The highest number of alleles was at 
MAF70 (34). At this locus, 23 alleles were detected in 
Ujmuqin sheep, and 22 alleles in other European sheep 



breeds (Dongyan et al, 2007; Lawson-Handley et al, 
2007). This shows the pronounced genetic difference 
among the flocks studied. For the loci BM1314, BM6506, 
MAF65, MCM527 and OarCP20, Cinkulov et al. (2003) 
found 1 1, 5, 5, 9 and 8 alleles, respectively, in the Serbian 
Cokanski Tsigai breed, whereas 5, 2, 4, 6 and 2 alleles, re- 
spectively, were observed at the same loci in our study. Ex- 
pected heterozygosity values over the populations ranged 
from 0.580 (MAF65) to 0.867 (MAF70), the latter compa- 
rable to the average heterozygosity value of 0.86 in Spanish 
sheep (Arranz et al. 2001). At the MAF65 locus, mean ob- 
served heterozygosity was 0.682 in Baltic sheep, 0.70 in 
Portuguese coarse-wool sheep and 0.754 in Alpine sheep 
breeds (Grigaliunaite et al, 2003; Dalvit et al, 2008; San- 
tos-Silva et al, 2008,). The mean number of alleles ranged 
from 6.3 (Rodopski Tsigai) to 8.6 (White Maritza), and was 
higher than that detected by Peter et al. (2007) in 57 Euro- 
pean sheep breeds, and Tapio et al. (2005) in northern Eu- 
ropean breeds. Santos-Silva et al. (2008) and Wafula et al. 
(2005) found similar values in Portuguese coarse-wool and 
West African sheep breeds (6.4-9.1, 6.5-7.4). The genetic 
diversity value (0,736) found in Bulgarian sheep is similar 
to that published by Grigaliunaite et al. (2003) for Baltic 
sheep (0.7 12), Oliveira et al. (2005) for Bordaleira de Entre 
Douro e Minho sheep (0.74) and Alvarez et al. (2004) for 
Latxa sheep (0.77). 

An analysis of Nei genetic distance between Bulgar- 
ian Tsigais, Maritzas and Pleven Black Heads indicated 
that the two Tsigais and the two Maritzas are closely re- 
lated, whereas Pleven Backheads are genetically unique. A 
genetic distance value of 0.05 implies moderate differentia- 
tion between two breeds (Hard, 1980). In our study, all ge- 
netic distance values for pairs of breeds were higher than 
0.05, which may indicate certain differences within their 
genetic structure. Considering the relatives of various 
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sheep breeds, Kukovics and Kume (2006) suggested that 
Pleven Blackheads might belong to the Ruda group. The 
two Bulgarian Tsigai breeds are completely different from 
other Tsigai variants, when considering phenotypic charac- 
teristics, as body measures, and face and leg coloring. 
These still carry the wool traits received from Soviet Tsigai 
breeds that were developed by crossing Tsigai sheep with 
British longwool breeds (Kukovics and Javor, 2002; Di- 
mov, 2006). 

Based on our results, we can affirm that Pleven 
Blackheads genetically differ from both the Bulgarian 
Tsigais and Maritza sheep, although they are more closely 
related to the latter. This provides a certain support to the 
hypothesis of Kukovics and Kume (2006), although the 
close relationship of Pleven Blackheads to Tsigai crossbred 
sheep (Maritzas) is also apparent. The black head and legs, 
as well as wool characteristics are similar to those of the 
original Tsigai sheep, but the large nose is visible evidence 
of an alternative original breed. 

In summary, the four breeds, two Bulgarian Tsigai 
and two Maritza, were found to be genetically closely re- 
lated to each other, although there was a significant dis- 
tance between the two groups. Pleven Blackhead sheep 
were found, as a whole, to be genetically distinct, although 
more closely related to the Maritzas. 
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